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ABSTRACT 


Intcrannual  and  climatological  variations  of  wind  stress  drag  coefficient  (Cn)  are  examined  over  the 
global  ocean  from  1908  to  2004.  Here  C „  is  calculated  using  high  temporal  resolution  (3-  and  6-hourly) 
surface  atmospheric  variables  from  two  datasets:  1 )  the  40-vr  European  Centre  for  Medium-Range  Weather 
Forecasts  (ECMWF)  Re-Analysis  (ERATO)  and  2)  the  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS).  The  stability-dependent  Cn  algorithm  applied  to  both  datasets  gives  almost  identical 
values  over  most  of  the  global  ocean,  confirming  the  validity  of  results.  Overall,  major  findings  of  this  paper 
are  as  follows:  1)  the  C n  value  can  change  significantly  (e.g.,  >50%)  on  12-hourly  lime  scales  around  the 
Kuroshio  and  Gulf  Stream  current  systems:  2)  there  is  strong  seasonal  variability  in  C,„  but  there  is  not 
much  intcrannual  change  in  the  spatial  variability  for  a  given  month;  3)  a  global  mean  Cn  **  1.25  x  10 
is  found  in  all  months,  while  Cn  ^  1.5  x  10  Ms  prevalent  over  the  North  Pacific  and  North  Atlantic  Oceans 
and  in  southern  high-latitude  regions  as  well,  and  Cn  ^  1.0  x  10  Ms  typical  in  the  eastern  equatorial  Pacific 
cold  tongue;  and  4)  including  the  effects  of  air-sea  stability  on  Cn  generally  causes  an  increase  of  >20%  in 
comparison  to  the  one  calculated  based  on  neutral  conditions  in  the  tropical  regions.  Finally,  spatially  and 
temporally  varying  Cn  fields  are  therefore  needed  for  a  variety  of  climate  and  air-sea  interaction  studies. 


I.  Introduction  and  motivation 

The  momentum  exchange  between  the  atmosphere 
and  ocean  through  wind  stress  is  of  importance  for 
many  purposes,  including  air-sea  interaction  studies, 
climate  studies,  ocean  modeling,  and  ocean  prediction. 
Wind  stress  is  typically  obtained  from  bulk  parameter- 
izations  that  estimate  turbulent  fluxes  using  standard 
meteorological  data  (e.g..  Fairall  el  al.  2003).  In  par¬ 
ticular,  the  total  wind  stress  magnitude  (r)  at  the  ocean 
surface  can  be  calculated  from  the  square  of  the  wind 
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speed  at  10  m  above  the  sea  surface  (F,,)*  the  density  of 
air  (p„),  and  a  dimensionless  drag  coefficient  (C,>)  with 
T  -  Pa  C„  Vi 

Previously,  attention  has  been  given  to  construction 
of  wind  stress  climatologies  over  the  global  ocean 
(Chelton  et  al.  1990)  and  experimental  analysis  of  Cn  at 
a  few  particular  locations  (Donelan  et  al.  1997).  Based 
on  the  authors'  knowledge  there  is  no  quantitative 
study  examining  the  spatial  and  temporal  variability  of 
Cn,  a  parameter  that  is  used  for  calculating  t.  over  the 
global  ocean.  Experimental  measurements  for  Cn  are 
rarely  available,  and  those  that  are  available  do  not 
have  sufficient  temporal  and  spatial  resolution  to  de¬ 
termine  its  global  distribution.  In  addition,  prior  to  the 
1990s  there  were  uncertainties  in  the  accuracy  of  near- 
surface  meteorological  variables  from  the  existing  ar¬ 
chived  numerical  weather  prediction  centers,  such  as 
the  European  Centre  for  Medium  Range  Weather 
Forecasts  (ECMWF)  and  the  National  Centers  for  En 
vironmenlal  Prediction  (NCF.P),  preventing  the  aeeu- 
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rate  calculation  of  CL)  based  on  bulk  variables.  For  ex¬ 
ample,  a  study  that  explores  climatological  mean  fea¬ 
tures  of  Cn  over  the  global  ocean  was  presented  by 
Trenberth  ct  al.  (1989).  That  study  was  limited  to  cal¬ 
culating  near-surface  stability  based  on  monthly  mean 
atmospheric  variables  (2.5°  x  2.5°  resolution)  at  the 
1000-mb  level  from  FCMWF  (1980-86),  which  was  the 
only  available  product  from  FCMWF  at  that  time.  That 
study  did  not  examine  diurnal  and  interannual  varia¬ 
tions  in  C n  either. 

The  quality  of  archived  products,  such  as  these  from 
ECMWF,  has  greatly  improved  since  that  time.  The 
amount  of  observational  data  has  greatly  increased  and 
better  assimilation  methods  have  been  developed. 
Thus,  it  is  now'  possible  to  better  determine  air-sea 
stability  and  to  investigate  spatial  and  temporal  varia¬ 
tions  (diurnal,  interannual,  and  climatological)  of  C,» 
over  the  global  ocean — the  major  goal  of  this  paper. 
This  will  be  accomplished  using  an  algorithm,  which 
fully  takes  the  air-sea  stability  into  account. 

A  study  that  examines  the  global  variability  of  Cn  is 
desirable  given  the  fact  that  the  use  of  an  inaccurate  Cn 
in  calculating  wind  stress  may  result  in  serious  errors. 
As  an  example,  ignoring  the  effects  of  water  vapor  flux 
in  the  parameterization  of  Cn  can  give  a  wind  stress 
value  that  is  ^6  times  less  than  its  actual  value  at  very 
low  wind  speeds  (Kara  et  al.  2005).  This  is  particularly 
true  on  short  time  scales  because  of  air-sea  stability 
(e.g.,  day  and  night).  In  other  w'ords,  averaging  stan¬ 
dard  meteorological  variables  over  a  day  may  generally 
give  neutral  conditions  (e.g.,  the  difference  between 
near-surface  air  and  sea  surface  temperature  is  almost 
zero  and  relative  humidity  is  100%).  This  is  generally 
not  the  case  for  daytime  versus  nighttime  conditions. 
For  this  reason,  C,,  needs  to  be  calculated  al  fine  tem¬ 
poral  resolution  (e.g.,  every  3  or  6  h)  and  then  averaged 
over  a  day. 

Given  the  need  for  proper  determination  of  stability- 
dependent  Cn  in  calculating  the  wind  stress  magnitude 
over  the  global  ocean,  the  main  focus  of  this  paper  is 
threefold:  1)  present  spatial  and  temporal  variations  in 
Cn  for  use  in  climate  studies,  2)  determine  regions 
where  there  is  strong/weak  seasonal  variability,  and  3) 
reveal  if  there  is  any  interannual  variability  in  Cn  or  if 
substituting  a  climatological  mean  may  be  appropriate 
over  the  different  regions  of  the  global  ocean. 

2.  Methods  and  data 

The  parameterization  of  Cn  is  still  an  active  field 
with  considerable  diversity  of  approaches  available  in 
the  literature  (e.g.,  Taylor  and  Yelland  2001;  Hwang 
2005).  Simply  slated,  the  drag  coefficient  is  usually  rep¬ 


resented  as  a  combination  of  a  wind  speed-dependent 
neutral  coefficient  (or,  equivalently,  a  roughness 
length),  a  hydrostatic  stability  dependence,  and  in  some 
cases  an  enhancement  at  low  wind  speeds  associated 
with  wind  gustiness  (e.g..  Fairall  et  al.  2003). 

Various  formulations  of  Cn  via  surface  roughness  are 
available.  They  are  based  on  friction  velocity  (rather 
than  V(l)  and  do  not  include  any  air-sea  stability  de¬ 
pendence  on  Cp  using  air-sea  temperature  difference 
or  vapor  mixing  ratio  values.  There  are  also  Cn  param¬ 
eterization  based  solely  on  Va  (e.g..  Trenberth  et  al. 
1989).  A  constant  CL)  value  has  been  used  in  many 
studies  (e.g.,  Kessler  and  Kleeman  2000;  Sura  et  al. 
2000;  Koracin  et  al.  2004).  Obviously,  w'ind  stresses  that 
are  calculated  using  such  constant  values  exclude  the 
significant  changes  in  magnitude  that  can  occur  due  to 
effects  of  air-sea  stability  in  C,>.  Several  formulations 
include  a  simple  form  for  stability  effects  on  Cn  based 
on  air-sea  temperature  difference  (e.g..  Smith  1988). 
However,  water  vapor  flux  is  also  an  important  param¬ 
eter  determining  stability,  especially  in  tropical  regions 
(Kara  et  al.  2005).  In  particular,  because  of  the  expo¬ 
nential  increase  of  saturation  vapor  pressure  with  tem¬ 
perature,  in  some  regions,  humidity  has  a  first-order 
effect  on  the  stability. 

For  the  reasons  mentioned  above,  a  parameteriza 
tion  that  takes  full  account  of  stability  in  calculating  Cn 
is  required.  The  detailed  parameterization  presented  in 
Kara  ct  al.  (2005)  is  used  here,  llie  formulation  is  based 
on  the  state-of-the-art  Coupled  Ocean-Atmosphere 
Response  Experiment  (COARE)  bulk  algorithm  (ver¬ 
sion  3.0).  which  employs  a  turbulence  theory  based  on 
the  iterative  estimations  of  the  scaling  variables  to  de¬ 
termine  stability-dependent  CD.  The  stability-depen¬ 
dent  C„,  as  used  in  this  paper,  is  expressed  as  simple 
polynomial  functions  of  air-sea  temperature  difference, 
using  air  temperature  at  10  m,  Va  at  10  m,  and  relative 
humidity  al  the  air-sea  interface  to  parameterize  air- 
sea  stability.  Because  of  deficiencies  in  the  COARF. 
algorithm  itself  at  high  winds  and  ongoing  debate,  Cn 
is  kept  constant  in  the  parameterization  for  winds 
>20  m  s-1. 

All  variables  for  calculating  Cn  over  the  global  ocean 
are  obtained  from  an  archived  gridded  (1  x  1”)  prod¬ 
uct — the  Fleet  Numerical  Meteorology  and  Oceanog¬ 
raphy  Center  (FNMOC)  Navy  Operational  Global  At¬ 
mospheric  Prediction  System  (NOGAPS;  Rosmond  et 
al.  2002).  NOGAPS  is  particularly  chosen  because  it 
provides  the  above-mentioned  atmospheric  variables  al 
high  temporal  frequency  (3-hourly),  a  critical  require¬ 
ment  to  take  air-sea  stability  into  account.  On  the  other 
hand,  in  section  4  we  will  also  use  similar  6-hourlv  at¬ 
mospheric  data  from  FCMWF  to  further  confirm  the 
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validity  of  Cf)  values  over  the  global  ocean.  The  drag 
CD  is  computed  from  1998  to  2004,  the  period  when 
suitable  archived  NOGAPS  data  are  available  as  of  this 
writing. 

Our  availability  of  3-hourly  surface  atmospheric  vari¬ 
ables  from  NOGAPS  started  in  2001;  thus  we  use 
6-hourly  data  from  1998  to  2000  and  3-hourly  data  from 
2001  to  2004.  The  Cn  is  first  calculated  at  each  3  (or  6)- 
hourly  time  interval  at  each  grid  point  over  the  global 
ocean.  Monthly  means  are  then  formed.  A  time  series 
of  1998-2004  may  not  be  long  enough  to  derive  a  cli¬ 
matological  mean  or  to  show  interannual  variability; 
however,  it  would  still  reveal  general  features  of  Cn 
over  the  global  ocean.  An  examination  of  Cn  per¬ 
formed  over  a  long  time  period  (e.g.,  1979-2002)  can  be 
accomplished  utilizing  atmospheric  variables  from  the 
ECMWF  Re-Analvsis,  a  topic  which  is  mentioned  in 
section  4. 


3.  Daily  cycle  of  stability-dependent  drag 

coefficient 

The  effect  of  stability  on  Cn  can  be  neglected  on 
longer  time  scales  (i.e.,  monthly).  This  is  justified  by  the 
fact  that  the  atmosphere  is  generally  neutral  or  slightly 
unstable  over  the  open  ocean.  Such  conditions  are  only 
true  for  wind  speed  values  between  6  and  25  ms  l,  as 
discussed  in  Bonekamp  et  al.  (2(K)2)  in  detail.  However, 
wind  speed  can  be  outside  this  range  over  many  regions 
of  the  global  ocean.  In  addition,  changes  in  stability 
through  air-sea  temperature  and  relative  humidity  dif¬ 
ferences  near  the  sea  surface  do  indeed  exist  when  con¬ 
sidering  the  diurnal  cycle  of  near-surface  atmospheric 
variables.  Phis  would  certainly  affect  C„,  a  stability- 
dependent  parameter. 

Thus,  a  few  particular  questions  arise  here:  how  vari¬ 
able  (spatially  and  temporally)  is  Cn  for  a  given  day? 
To  answer  these  questions  Cn  values  calculated  from 
the  formulation  of  Kara  et  al.  (2005)  using  data  from 
NOGAPS  (see  section  2)  are  analyzed  on  6  January 
2004.  I  his  day  is  chosen  just  for  illustrative  purposes. 
Figure  la  clearly  demonstrates  large  spatial  variability 
in  Cn  over  the  global  ocean. 

To  examine  the  temporal  variability  of  Cn  a  few 
zoom  regions,  surrounding  two  major  oceanic  current 
systems  (Gulf  Stream  and  Kuroshio)  and  the  western 
equatorial  Pacific  warm  pool,  are  selected.  The  Cn  is 
shown  at  00007  and  12007  (Fig.  lb).  The  increase  or 
decrease  in  Cn  from  00007  to  12007  can  be  larger  than 
50%  near  regions  around  the  Gulf  Stream  and  Kuro¬ 
shio  as  evident  from  the  ratio  values.  As  expected,  this 
variability  is  due  to  the  wind  speed  and  air-sea  stability 


changing  from  00007  to  12007.  For  example,  Vit  is  9.1 
(11.9  ms-1),  air-sea  temperature  difference  is  2.1 
(  5.3°C),  and  relative  humidity  is  89%  (66%)  at  38°N, 
70°W  in  the  Gulf  Stream  region  at  00007  ( 12007).  The 
corresponding  CD  X  10*  value  is  1.14  (1.56).  The  C',>  in 
the  western  equatorial  Pacific  has  much  less  variability 
than  that  in  the  Gulf  Stream  and  Kuroshio  current  sys¬ 
tems. 

Time  periods  (00007  and  12007)  were  chosen  ran¬ 
domly  in  order  to  investigate  changes  over  a  12-h  time 
period.  The  ratio  of  Cn  [i.e.,  C’/)(at  (K)()0Z)/C’/)(at 
12007)]  is  roughly  out  of  phase  in  both  the  Gulf  Stream 
and  Kuroshio  regions  if  one  considers  the  day/night 
cycle  (i.e.,  being  dav/night  in  the  Gulf  Stream  and  night/ 
day  in  the  Kuroshio).  To  capture  the  full  diurnal 
change,  the  comparison  in  each  region  is  also  made 
between  the  Cf)  map  at  the  3-hourly  sampling  time  clos¬ 
est  to  high  temperature  and  the  map  al  the  time  closest 
to  low  temperature.  The  resulting  ratio  values  clearly 
reveal  that  there  are  indeed  large  changes  in  Cn  from 
day  to  night  in  both  regions  (not  shown),  and  the  re¬ 
sulting  ratios  for  Cn  generally  resemble  those  shown  in 
Fig.  lb. 

Changes  in  Cn  mentioned  in  the  preceding  analysis 
were  for  a  particular  day.  Thus,  a  lot  of  the  differences 
could  be  temporal  but  not  necessarily  diurnal.  Thus,  we 
also  investigated  the  mean  diurnal  variability  of  Cn  by 
averaging  values  at  each  3-hourly  interval  from  00007 
to  12(X)Z  and  from  1200Z  to  2400Z  over  a  month.  The 
ratio  of  Cn  during  the  two  time  periods  is  generally 
found  to  be  within  ±5%  (±10%)  during  the  Northern 
Hemisphere  winter  (summer). 

Regarding  the  preceding  analysis,  we  need  to  empha¬ 
size  that  the  current  parameterization  of  Cn  ignored 
the  effects  of  ocean  currents  and  waves.  Such  factors 
are  particularly  important,  given  that  this  paper  focuses 
on  several  regions  experiencing  relatively  strong  cur¬ 
rents,  such  as  the  Gulf  Stream  and  the  Kuroshio.  For 
example,  it  is  well  known  that  especially  in  regions  of 
strong  currents,  it  is  not  simply  the  wind  speed  that  is 
important  but  the  vector  difference  in  wind  speed  and 
ocean  current  speed  and  waves.  In  addition,  possible 
impacts  of  ocean  currents  and  wand  waves  on  Cn  have 
been  discussed  in  various  studies  (e.g.,  Wuest  and 
Lorke  2003).  On  the  other  hand,  a  previous  study  based 
on  a  high  temporal  resolution  global  dataset  demon¬ 
strated  that  strong  ocean  currents  near  the  western 
boundaries  (Kuroshio  and  Gulf  Stream)  do  not  sub¬ 
stantially  influence  Cn  (Kara  et  al.  2007).  This  is  due  to 
the  fact  that  the  winds  and  currents  are  generally  not 
aligned  or  locally  correlated.  In  particular,  that  study 
revealed  that  the  combined  outcome  of  ocean  currents 
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(a)  Drag  coefficient  (Co  X  10’*)  over  the  global  ocean  on  G  Jan  2004  (00Z) 


(b)  Drag  coefficient  (Co  X  10’*)  at  three  zoom  regions  on  6  Jan  2004 

(00Z)  (12Z)  Ratio 


Fk..  1.  Values  of  drag  coefficient  (a)  for  the  global  ocean  at  0000Z  and  (b)  in  3  zoom  regions  numbered  as  1.  2,  and  3  (near  the 
Kuroshio.  Gulf  Stream,  and  equatorial  regions,  respectively)  at  0000Z  and  1200Z.  The  ratio  of  drag  coefficient  [C„(al  0000Z)/C/>(at 
1 2<X)Z )]  is  also  provided. 


and  waves  is  to  reduce  daily  Cn  by  <5%  within  those 
current  systems. 

4.  Intcrannual  variability  of  drag  coefficient 

Monthly  mean  C„  fields,  based  on  the  method  and 
data  described  (section  2),  are  shown  in  Fig.  2  over  the 


global  ocean  from  1998  to  2004.  For  simplicity,  we 
present  mean  Cn  fields  in  February,  August,  and  No¬ 
vember  for  each  year.  It  is  clear  that  there  is  no  noise 
nor  bull's-eyes  in  the  fields,  although  no  smoothing  was 
applied  to  the  Cn  fields.  Note  that  in  all  panels,  the 
regions  where  ice  exists  are  masked  out  since  Cn  over 


1998  to  2004  200  4  2003  2002  2001  2000  1999  1998 
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(ci )  February  (b)  August  (c)  November 


Cn  X  1 0'' 


Fk;.  2.  Monthly  mean  drag  coefficients  over  the  global  ocean  in  (a)  February,  (b)  August,  and  (c)  November.  They 
are  shown  from  1998  to  2004.  along  with  mean  drag  coefficients  calculated  over  the  7-yr  lime  period  (bottom  row  ). 
All  values  arc  in  l()\  Regions  where  ice  is  present  are  shown  in  gray. 
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ice  is  quite  different  from  that  over  the  sea.  and  it  is  not 
the  purpose  of  this  paper  to  discuss  sea  ice  drag  coef¬ 
ficients.  The  ice-free  regions  over  the  global  ocean  (ice 
covered  shown  in  gray)  are  determined  using  an  ice- 
land  mask  based  on  data  described  in  Reynolds  et  al. 
(2002). 

The  most  significant  feature  evident  from  fields  is 
that  Cn  changes  month  by  month  during  1998-2004. 
This  suggests  that  spatial  and  seasonal  variations  in  Cn 
need  to  be  taken  into  account  in  climate  applications, 
especially  in  the  high  northern  and  southern  latitudes. 
On  the  contrary,  there  is  little  interannual  variability 
for  a  given  month.  The  Cn  values  that  are  as  large  as 
1.5  x  10”1  in  February  become  as  low  as  *=<1.1  x 
10  '  in  August  over  most  of  the  North  Pacific  and 
North  Atlantic  Oceans,  including  the  regions  around 
the  Kuroshio  and  Gulf  Stream  current  systems.  In  fact, 
a  close  examination  of  spatial  fields  near  the  Kuroshio 
region  reveals  that  there  is  an  **>50%  change  or  more 
in  Cn  from  February  (si. 6  x  10  *)  to  August  (^1.1  x 
10  3).  The  Cp  value  increases  again  in  November.  This 
is  not  surprising  given  the  fact  that  there  are  even  larger 
differences  for  a  given  day  (see  Fig.  1).  The  strong  sea¬ 
sonality  of  Cp  is  also  evident  at  high  southern  latitudes. 

Previously,  Trenberth  et  al.  (1989)  also  noted  similar 
results,  in  that  Cp  at  high  latitudes  is  generally  larger 
than  in  the  Tropics.  However.  Cn  values  as  low  as  <1  X 
10  \  which  are  seen  in  the  eastern  equatorial  Pacific 
cold  tongue  in  August  and  partially  in  November  (Fig. 
2),  were  not  detected  in  their  study.  This  is  due  mainly 
to  the  different  gridded  data  and  resolution  used  for 
calculating  Cn  [the  coarse-resolution  2.5°  X  2.5° 
ECMWF  1000-mb  data  (1980-86)  versus  the  relatively 
fine  1.0°  X  1.0°  NOGAPS  surface  data  (1998-2004) 
used  in  this  paper].  The  winds  at  1000  mb  can  also  be 
very  different  from  the  winds  10  m  above  the  sea  sur¬ 
face.  However,  those  were  the  only  data  available  at 
that  time.  In  addition,  there  were  known  deficiencies  in 
the  earlier  ECMWF  data,  especially  in  tropical  regions, 
as  indicated  in  their  study. 

One  important  remark  here  addresses  an  aspect  of 
uncertainty  in  the  value  of  Cn  presented  in  Fig.  2.  As 
mentioned  previously,  Cn  was  calculated  using  atmo¬ 
spheric  variables  (air  temperature  and  wind  speed  at  10 
m,  sea  surface  temperature,  and  relative  humidity  at  the 
air-sea  interface)  from  NOGAPS.  Possible  inaccura¬ 
cies  that  may  exist  in  these  variables  can  be  associated 
with  errors  in  Cn.  Thus,  one  might  argue  that  there  is 
not  nearly  enough  information  to  make  an  estimate  as 
to  whether  the  errors  in  atmospheric  variables  from 
NOGAPS  would  outweigh  any  signal  in  Cn  variations 
discussed  in  this  paper. 

For  the  reasons  mentioned  above,  we  also  calculate 


Cp  by  substituting  the  similar  atmospheric  variables 
from  another  global  dataset,  the  40-yr  ECMWF  Re- 
Analysis  (ERA-40)  from  the  years  1979  to  2002  (Kall- 
berg  et  al.  2004).  For  simplicity.  Cn  values  calculated 
from  NOGAPS  and  ECMWF  data  are  compared  for 
the  months  of  January  and  June  (Fig.  3).  This  is  done 
for  two  different  years,  2000  and  2001 .  to  investigate  the 
consistency  of  the  Cn  comparison  in  different  years. 
Clearly,  there  are  no  significant  differences  between  Cn 
values  calculated  from  NOGAPS  and  ERA-40,  as  evi¬ 
denced  by  Cp  ratios  close  to  1  over  most  of  the  global 
ocean  (Fig.  3).  The  seasonal  variability  of  Cn  is  prop¬ 
erly  represented  using  either  one  of  the  datasets.  Com¬ 
parisons  of  atmospheric  variables  (e.g.,  wind  speed,  air 
temperature,  etc.)  between  NOGAPS  and  ERA-40  arc 
beyond  the  scope  of  this  paper.  However,  as  demon¬ 
strated  here,  both  datasets  give  almost  identical  Cn  val¬ 
ues  over  most  of  the  global  ocean,  implying  close  agree¬ 
ment  in  surface  atmospheric  variables  between  the  two. 

Finally,  we  also  investigate  Ihe  impact  of  stability  on 
Cp.  It  is  possible  that  a  user  would  like  to  compute  Cn 
based  on  only  wind  speed  because  air-sea  temperatures 
and  humidity  may  not  be  available  or  missing  in  an  in 
situ  dataset.  Thus,  the  user  would  like  to  know'  the  ac¬ 
curacy  penalty  when  only  wind  speed  is  used  for  calcu¬ 
lating  Cp.  To  answer  this  question,  more  calculations 
are  performed.  Essentially,  we  use  6-hourly  atmo¬ 
spheric  variables  from  ERA-40,  spanning  1979  through 
2002.  We  calculated  Cn  using  (i)  the  air-sea  tempera¬ 
ture  difference  of  zero  along  with  relative  humidity  of 
100  (i.e..  neutral  conditions)  and  (ii)  the  actual  air-sea 
temperature  difference  and  relative  humidity  (full  sta¬ 
bility).  Both  (i)  and  (ii)  are  done  at  each  grid  over  the 
global  ocean. 

We  first  calculate  Cn  and  wind  stress  at  each  6-hourly 
time  interval  and  then  form  monthly  mean  Cn  and  w'ind 
stress  for  each  year.  Finally,  a  monthly  mean  climatol¬ 
ogy  is  formed  for  both  quantities  over  1979-2002.  This 
is  done  to  provide  an  easier  presentation  of  the  results. 
The  difference  (long-term  mean  bias)  between  Cn  cal¬ 
culated  based  on  full  stability  (ii)  and  that  calculated 
based  on  neutral  conditions  (i)  clearly  demonstrates  a 
significant  effect  of  air-sea  stability  on  Cn  even  on 
monthly  time  scales  (Fig.  4a).  The  Cn  including  effects 
of  the  stability  is  generally  larger  over  most  of  the  glo¬ 
bal  ocean.  If  one  computed  Cn  based  on  neutral  stabil¬ 
ity  as  opposed  to  the  full  stability,  errors  could  be  very 
large  and  might  be  larger  than  20%  especially  in  the 
tropical  regions  (Fig.  4b). 

5.  Summary  and  conclusions 

There  are  no  detailed  studies  examining  the  variabil¬ 
ity  of  C n  on  shorter  (e.g.,  diurnal  and  daily)  and  longer 
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Co  X  10'*  (NOGAPS)  CD  x  10'*  (ERA  40)  Ratio 


Fio.  3.  Values  of  drag  coefficient  calculated  using  surface  atmospheric  variables  from  NOGAPS  and  ERA-40  as  explained  in  the  text. 
ITic  ratio  of  drag  coefficient  [Cr>(NOGAPS)/Cn(ERA-40)]  is  also  provided  and  white  represents  ratios  between  0.95  and  1.05. 


(monthly,  intcrannual,  and  climatological)  time  scales 
over  the  global  ocean,  even  though  there  may  exist  field 
studies  at  specific  time  periods  and  certain  places.  ITiis 
information  about  spatially  and  temporally  varying  CD 
over  the  global  ocean  is  therefore  long  overdue  in  the 
literature.  The  aim  of  this  paper  is  to  fill  in  this  gap. 
Such  information  is  greatly  desired  by  not  only  climate 
modelers  but  also  modelers  dealing  with  coupled  atmo¬ 
sphere-ocean  and  ocean-only  studies.  In  many  cases, 
researchers  dealing  with  such  studies  would  like  to  see 
maps  of  these  fields,  as  provided  in  this  paper,  to  un¬ 
derstand  the  variability  of  this  particular  quantity,  es¬ 
pecially  for  calculating  wind  stress. 

In  this  paper,  we  use  high  temporal  (e.g..  3-hourly) 
and  fine  spatial  (1.0°  x  1.0°)  resolution  near-surface 
atmospheric  variables  from  NOG  AES  to  determine  sta¬ 
bility-dependent  Cf)  on  both  interannual  and  climato¬ 
logical  time  scales.  An  examination  of  Ciy  fields  around 
the  major  oceanic  current  systems  (e.g.,  the  Kuroshio 


and  the  Gulf  Stream)  demonstrates  spatial  variations 
not  only  for  day  versus  night  conditions  but  also  for 
monthly  time  scales  from  1998  to  2004.  For  example,  a 
CL)  value  of  1.0  X  10  '  near  the  Gulf  Stream  at  0000Z 
can  become  as  large  as  1.8  x  10  *  at  1200/.,  implying  a 
diurnal  change  *^80%.  This  indicates  the  importance  of 
including  stability  effects  in  calculating  Ct),  and  thereby 
wind  stress  since  air-sea  temperature  differences,  rela¬ 
tive  humidity  in  addition  to  wind  speed  can  be  quite 
variable  at  those  time  periods.  Outside  the  tropical  re¬ 
gions,  CD  values  are  quite  different  from  ***1.2  X  10  ' 
and  can  even  be  ^1.6  x  10  '  with  strong  seasonal  vari¬ 
ability — an  important  feature  of  Cn  that  should  be 
taken  into  account  in  air-sea  interaction  and  climate 
modeling  studies. 

We  also  calculated  monthly  mean  Cn  fields  using 
6-hourly  atmospheric  variables  from  the  1.125°  x 
1.125°  ERA-40  product.  The  Cn  values  were  found  to 
be  similar  to  those  calculated  from  NOGAPS.  indicat- 
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(a)  Mean  bias  in  “neutral44  drag  coefficient  (Co  X  10,t)  (1979  2002) 
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Fk..  4.  (a)  Difference  between  ihe  drag  coefficient  calculated  using  full  air-sea  stability  dependence 
and  the  one  calculated  using  neutral  conditions.  The  latter  is  subtracted  from  the  former.  As  explained 
in  the  text,  differences  in  C  n  arc  calculated  at  each  grid  point  over  the  global  ocean  during  1979  2002. 
(b)  Percentage  ratio  of  the  drag  coefficient  computed  using  full  air  sea  stability  to  the  one  computed 
using  neutral  conditions  over  the  global  ocean  during  1979-2002. 


ing  l he  robustness  of  the  fields  presented  in  this  paper. 
One  important  conclusion  arising  from  this  paper  is 
that  ignoring  the  effects  of  air-sea  stability  in  comput¬ 
ing  Cp  generally  results  in  >10%  error  in  comparison 


to  the  actual  error  that  lakes  full  stability  into  the  ac¬ 
count.  Such  errors  are  even  larger  (typically  >20%)  in 
the  tropical  regions.  Finally,  an  examination  of  Cn  over 
the  global  ocean  using  other  parameterizations  may 
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yield  slightly  different  results,  but  sueh  an  investigation 
is  beyond  the  seope  of  this  paper. 
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